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ABSTRACT 

With the Spzfaer-IRS^ we have obtained the first mid-IR spectroscopy of NGC 7129, in the unusuaUy 
strong outflow and in a ridge of II2 emission near the Herbig Be star LkHa 234. The UV radiation 
field strength is estimated from PAH band intensities in the H2 ridge, and found to be comparable to 
that of NGC 7023. From the rotational H2 emission lines we have deduced aperture average excitation 
temperatures and column densities in the two regions, finding the H2 ridge values to be consistent 
with pumping by UV fluorescence, but also comparable to warm-gas regions of Cep A that form H2 in 
non-dissociative C-shocks. The H2 emission in the outflow is consistent with formation by collisional 
excitation in J-shocks, with shock velocities of 10 - 30 km sec~^. A photodissociating component 
may be present in the outflow, by similarity of S(0) line intensities in both regions. There is no 
indication of warm dust in the outflow. We also present the flrst 16 /xm imaging of a Galactic nebula 
using the unique imaging capabilities of the IRS, and combine with ground-based 2.12 /^m (H2 1 - 
S(l)) imaging. Candidate pre-main sequence objects are clearly evident in these data. We also 
find extended emission not previously observed around the young B star BD -|-65°1638, '^ 22" across, 
showing that the region is not free of material as otherwise inferred by recent high angular resolution 
mapping at submillimeter wavelengths. The presence of this material complicates interpretation of 
the surrounding CO cavity and origin(s) of the photodissociated region, and further spectroscopic 
observations are needed to characterize its nature. 

Subject headings: infrared: ISM — ISM: jets and outflows — ISM: individual (NGC 7129) — reflection 
nebulae — stars: formation — stars: pre-main sequence — infrared: stars 



1. INTRODUCTION 

NGC 7129 in Cepheus is a prototype young star- 
forming region and well-known reflection nebula, domi- 
nated by two young intermediate-mass [M ~ 3 — lOM©) 
B stars, and the unusual Herbig Be star LkHa 234. The 
winds from the early-type B stars, BD +65°1637 and 
BD -|-65°1638 (B2-B4; Strom et al.l972; Hillenbrand et 
al. 1992), but chiefly from the latter which is proba- 
bly the older (Racine 1968), have created a cavity of 
CO emission near the center of the nebula, and dense 
ridges of molecular material to the east, south, and west 
(MitcheU & Matthews 1994; Fuente et al. 1998; Eisloffel 
2000; Miskolczki et al. 2001). Atomic and molecular 
lines observed in an 80" beam with the Infrared Space 
Observatory Long Wavelength Spectrometer (LWS) in- 
dicate that the emission originates from photodominated 
excitations, possibly in high density or clumpy environ- 
ments as well (Lorenzetti et al. 1999; Giannini et al. 
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1999). LkHa 234 is estimated to be the youngest of 
these dominant stellar sources (Hillenbrand et al. 1992), 
probably formed in response to the compression of mate- 
rial on the northern edge (Mitchell & Matthews), and is 
in a center of star forming activity that includes several 
pre-main sequence (PMS) objects, observed at millimeter 
wavelengths with high angular resolution (Fuente et al. 
2001). A low velocity (~10 km sec~^) outflow extends 
over ~ 12' from LkHa 234 (Edwards & SneU 1983), pos- 
sibly as far as ~ 22' or some 8 pc as projected onto the 
sky (McCroarty, Ray, & Bally 2004) at a distance of 1.25 
kpc to NGC 7129 (Shevchenko & Yakubov 1989), and an 
optical jet extends into the cavity (Ray et al. 1990). The 
outflow and jet likely originate from a PMS object offset 
from LkHa 234, as suggested by Cabrit et al. (1997) and 
then localized by Fuente et al. (2001) to a source ~ 3". 5 
to the NW. 

In this paper we present the first mid-infrared spectra 
of two regions, one in the outflow, and one in the ridge 
of H2 emission near LkHa 234, with the aim of exploring 
the excitation conditions in the environment surround- 
ing this young star in its early stages of evolution. We 
estimate the intensity of the UV radiation field from the 
strengths of the polyaromatic hydrocarbon (PAH) band 
strengths, and determine the excitation mechanism re- 
sponsible for H2 line emission. Imaging the nebula at 
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Fig. 1. — Position of the IRS spectral slits on NGC 7129, 11 November 2003. The slits are Short Low (white or black), Long Low 
(yellow), Short High (green), and Long High (blue). The 14' X 12' image is an IRAC mosaic with the 3.6 (blue), 4.5 (green), 5.8 (orange), 
and 8.0 /im (red) filters (see Gutcrmuth ot al. 2004). The blue dotted lino indicates where IRS 16/^m and ground-based 2.12 fim mapping 
have been obtained (Fig.l^. 



2.12 /Ltm (H2 1-0 S(l)) and 16 ^m (continuum plus H2 
- S(l)) gives us a larger scale view on the spatial 
distribution of molecular gas. 

2. SPECTROSCOPY AND IMAGING 

2.1. IRS Spectroscopy 

NGC 7129 was observed on 11 November 2003 with 
the Infrared Spectrograph (IRS, described by Houck et 
al. 2004, this volume) onboard the Spitzer Space Tele- 
scope (Werner et al. 2004a, this volume). A spectral map 
was performed with the slits centered on the II2 ridge 
in the photodissocated region (PDR), at 2lh43'"04.82" 
-H66°Q6'27".7 (all positions J2000.0), 29". 5 arcsec SW of 
LkHa 234. The map also includes a position in the strong 
molecular outflow, 2lh43'"19.8'' -K66°07'50".7. Our anal- 
yses are concentrated on these two positions, nearest to 
LkHa 234, and where the fullest wavelength coverage 
among the separate slits is most spatially confined (as 
described further below) . All slit positions on the nebula 



are shown in Fig.Q using the Early Release Observation^ 
of the nebula, obtained with the Infrared Array Camera 
(IRAC; Fazio et al. 2004, this issue) at 3.6, 4.5, 5.8, and 
8.0 /im. 

Spectroscopy was obtained with each of the IRS mod- 
ules, using total exposure times of 48 seconds with the 
Short Low module (SL, spanning 5.3 - 14.5 /xm at spec- 
tral resolution R = A/AA ~ 80 - 120), 90 seconds with 
Short High (SH, 10 - 19.7 /um, R ~ 700), and 120 seconds 
with Long High (LH, 18.8 - 37 /um, R ~ 700) at each slit 
position. Long Low (LL, 14 - 38 /im, i? ~ 80 - 120) 
data were also obtained, providing an internal check of 
the consistency with overlapping SH and LH data. The 
relative spectrophotometry is consistent to within 10% 
over the same sky area, and the LL data contain no new 
spectral information, rather degraded at the lower spec- 

^ http://www.spitzer.caltech.edu/Media/releases/ssc2004- 

02/ssc2004-02a.shtml. See also Gutermuth et al. (2004) and 
MuzeroUe et al. (2004). 
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tral resolution; therefore we present only the SH and LH 
data here. The raw data have been processed in the SSC 
S9.1 pipeline to Basic Calibrated Data (BCD) products, 
and then extracted to 1-D spectra using the offline ver- 
sion of the post-BCD pipeline. A region of 3".6x 13". 5 on 
the central map position was extracted from the SL data 
at 10 /zm. The full SH and LH slits were extracted and 
then scaled to the SL spectrum. Two or three individual 
spectra over each wavelength range were obtained at each 
extracted position, and these were flux-calibrated, coad- 
ded, and sigma-clipped to produce a single spectrum of 
the H2 ridge spanning 5.3-37 /zm, and an outflow spec- 
trum over 7.5 - 37 /im. Absolute fluxes are estimated to 
be uncertain by ^20% overall (Decin et al. 2004). 

2.2. Ground-Based 2.12 ^m and IRS 16 ^m Imaging 

The H2 1 - S(l) 2.l2fim data were obtained on 2003 
November 10 - 12 with the 3.6m telescope at the Canada- 
France-Hawaii Telescope (CFHT) using the CFHTIR 
camera. The camera provides a 3'. 6 x 3'. 6 field of view 
and a scale of 0".211 pixel" ^. The data were dark sub- 
tracted, flatfielded, and sky corrected in IRAF, using the 
special near-IR data reduction package DIMSUM"^ . The 
H2 image has a total integration time of 18 minutes, and 
been resampled to match the pixel scale of the 16^m IRS 
image. 

Mid-IR imaging of NGC 7129 was obtained with the 
IRS 16 /im camera on 10 October 2003 and 29 Novem- 
ber 2003. Individual data frames were reduced in the 
S9.1 SSC pipehne to BCD products, and then registered 
based on reconstructed pointing, which is accurate to 
~0.7" (1-cr radial). The 16 /im and 2.12 /im mosaic is 
shown in Fig. 13 Several pixels are saturated at 16 /im 
by the three identified objects. In this field, MuzeroUe 
et al. (2004, this issue) have proposed several PMS ob- 
jects on the basis of the colors extracted from IRAC and 
MIPS, three of which are clearly evident by their 16 /tm 
emission in Fig. 13 White arrows identify candidate Class 
O/I objects, the yellow arrow corresponds to a Class II 
object. At least one other PMS object in this field is 
proposed by MuzeroUe et al., but is not evident by a 16 
/im thermal excess. 

Fig. 121 also shows an arc of 16 /im emission in the vicin- 
ity of BD -1-65° 1638, almost concentric with the eastern 
H2 ridge. This structure is not discerned in IRAC or 
MIPS imaging (cf. Gutermuth et al. 2004 and MuzeroUe 
et al. 2004), possibly as a consequence of the tempera- 
ture and surface brightness of the material. At moderate 
angular resolution (1'), Matthews et al. (2003) iden- 
tify a dense concentration of neutral hydrogen (21 cm) 
coincident with BD -|-65°1638, as well as an extensive 
region of photodissociated H2 accompanied by a small 
H II region around this star. Miskolczi et al. (2001) have 
proposed that the cavity of molecular material around 
BD +65° 1638 has been created by photodissociation of 
H2 into neutral hydrogen, and that the gas and radia- 
tion pressure from within the cavity are responsible for 
sweeping material into a shell that appears as the H2 
ridge. The presence of this concentrated structure (ac- 
knowledging some saturation of an area of radius ^ 2". 5 

^ The Deep Infrared Mosaicing Software package of IRAF 
scripts, developed by P. Eisenhardt, M. Dickinson, A. Stanford, 
and J. Ward, available at ftp://iraf.noao.edu/contrib/dimsuniV2J 



centered on BD +65° 1638) clearly shows that material, 
either molecular gas or dust, in not entirely evacuated 
from this region. Further spectroscopy is needed deter- 
mine its properties. 



TABLE I 
Rotational H2 - observed fluxes, column densities, and 

EXCITATION TEMPERATURES^. 
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8.75 (0.77) 
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S{2) 12.28 


1.20 (0.03) 


40.84 
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0.54 (0.02) 


40.04 


1682 


S{I) 17.03 


1.70 (0.04) 


42.41 


844 


0.54 (0.02) 


41.27 


844 


S(0) 28.22 


0.37 (0.02) 


45.60 


510 


0.33 (0.01) 


45.49 


510 






A'tot 
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A^tot 
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1.17 (0.04) 


658 (9) 




13.4 (0.7) 


248 (12) 



^Fluxes in units of lO"'^'^ W cm~^ sr" 
parentheses. Excitation temperature Te- 
em' 



uncertainties are in 
is in K, and total col- 



umn density density A^tot is in lO"'® cm ^. 
''Severely blended with neighboring PAH band. 
"^Blended with neighboring PAH band. 

3. SPECTROSCOPIC RESULTS 

The IRS spectra of the PDR/H2 ridge and outflow are 
presented in Fig. El as labeled. The spectra are not off- 
set from each other. Sky background levels indicated by 
the triangles at 10, 24, and 30 /im are estimated from a 
Spitzer-centvic, COBE/DIRBE-based^ model of the ther- 
mal emission by zodiacal and ISM dust in the line of sight 
to NGC 7129, on its date of observation. At these levels, 
we surmise that the outflow contains negligible amounts 
of warm dust. 

The PDR spectrum exhibits the usual set of PAH 
bands, at peak positions of 6.22, 7.63, 8.63, 11.22, 12.75, 
and 16.46 /tm, dominated by the 7.63 and 11.22 /im 
bands. Secondary peaks that may be associated with 
PAH emission are present at 11.00, 11.99, 13.55, 14.24, 
and possibly 17.37 /tm. These features are not previously 
identified, to our knowledge, and Werner et al. (2004b, 
this volume) identify other new subsidiary bands in IRS 
spectra of the NGC 7023 reflection nebula, including the 
17.4 /tm feature. 

The strength of the PAH band emission allows us to 
estimate the intensity of the UV radiation fleld in the 
vicinity of LkHa 234. Using a two-component blackbody 
with temperatures 75 K and 320 K and unity emissiv- 
ity to approximate the thermal continuum, we measure 
strengths of the 6.22, 7.63, 11.22, and 12.75 /im bands in 
the range of 6.0 — 30 x 10~^ W m^^ sr^^. The observa- 
tional trend of PAH band intensities plotted against UV 
radiation field density by Boulanger et al. (1998) for a 
set of PDRs implies a field strength of 800 to 900 times 
the solar value. This is close to the value of NGC 7023, 
which is measured by Boulanger et al. from a region 
of maximum emission in the PDR, in proximity to the 
PMS Hcrbig B2-B3 star HD 200775 which is similar to 
our PDR location near BD+65° 1638, without consider- 
ation in either case for shielding by circumstellar mate- 
rial. And like NGC 7023, the PAH bands are lower in 

* Avaialable in the SPOT observations planning tool. 
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Fig. 2. — IRS 16 /im (red) and CFHT 2.12 /xm (green) mosaic of NGC 7129. The bandpass of the 16/.im camera covers 13.3 - 18.7 fira. 
Individual frames cover 54" X 82" with a pixel scale of 1".8. A mosaic of 25 positions, with 140 seconds of integration time per pixel on the 
sky, was obtained to cover an area of approximately 4'. 3 X 2'. 5. Regions of spectroscopic analysis are indicated by the circles with dashed 
lines. Candidate PMS objects are indicated by the yellow and white arrows. A region of ~ 80" X 85" around BD +65° 1638 at 16/im is 
shown in the inset, to bring out details of structured emission around the star, and the region of relative maximum surface brightness in 
the H2 ridge. 



intensity at this UV field strength than given by a sim- 
ple power law extrapolation from PDRs with lower field 
strengths. This could be due to a lower abundance of 
small dust grains in high radiation fields, as suggested 
by Boulanger et al. 

The FDR and outflow spectra both exhibit pure rota- 
tional lines of H2 , arising from optically thin quadrupole 
transitions that are readily thcrmalized at moderate vol- 
ume densities. This allows us to estimate H2 column 
densities and excitation temperatures under the assump- 
tion of LTE (e.g., Gredel 1994), in order to gain in- 
sight into the excitation conditions of the line forming 
regions. We follow excitation diagram methods devel- 
oped by Burton (1992) and Gredel (1994), and applied 
to comparable environments (as in, e.g., Cepheus E by 
Moro-Martin et al. 2001). The column density for tran- 
sition i is Nj = AirFj/nEjgjAj, where Ej is the en- 
ergy of the upper level, gj is the statistical weight of 
ortho- and para-transitions, and Aj is the transition 
probability. The line column density is related to the 
total column density and excitation temperature Tox by 
\n{Nj/gj) = ln{Ntot/Q) — Ej/{kTc^), and may be solved 
for with a Boltzmann excitation diagram. The partition 
function Q was determined with an ortho-para ratio of 
3 under the assumption of thermal equilibrium (Moro- 
Martin et al. 2001). For the observed line fluxes per 
unit solid angle Fj/Q, we we have extracted over regions 
of 40 arcsec^ and 62 arcsec^ in the 5.3 - 14 /im and 
14 - 35 fim ranges, respectively. In these areas we must 



view our measurements as spatially integrated. The com- 
puted values for ln(Nj/gj), Eup/K along with the final 
estimates of the total H2 column density and excitation 
temperature for the FDR and outflow are listed in Table 
1. The observed line fluxes per unit solid angle Fj/Q, 
1-cr measurement uncertainties, and computed values of 
\n{Nj/gj), Eup/k along with final estimates of the to- 
tal H2 column density and excitation temperature for 
the FDR and outflow are listed in Table ^ The lines 
were measured with the line- fitting utilities in the ISAP^, 
which allows multiple fits of Lorentzian and Gaussian 
profiles for deblending, suited to the PAH features (e.g., 
Boulanger et al. 1998) and the unresolved H2 lines. This 
permits the H2 S(4) and S(5) fluxes to be separated from 
the 7.6 and 8.6 /im FAH bands. The FDR line intensi- 
ties are corrected for reddening using the Draine & Lee 
(1984) extinction law, and Ay = 2.6 mag estimated by 
Hamann & Persson (1992) towards LkHa 234. The ex- 
tinction across the nebula probably varies widely, and 
the ISM component is likely to be much lower along the 
line of sight to our II2 ridge observation, but the differ- 
ences to the S(0) - S(2) lines will be small (< 5%), and 
measurement uncertainties of the S(3) - S(5) lines arise 
predominantly from the blending with the PAH features 
and the flux calibration. The S(6)/S(5) or S(4)/S(5) Hne 
ratios could be also affected by an ortho-para ratio differ- 

^ ISO Spectral Analysis Package: 

|http://www.ipac.caltech.edu/iso/isap/isap.html| 
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Fig. 3. — IRS spectra of the outflow and PDR near LkHa 234, as indicated. Data are on their natural flux scales, without offsetting. 
Diamonds represent emission levels of the sky (ISM and zodiacal light) background. The atomic nebular and H2 rotational are identifled, 
where dashed lines indicate where H2 lines in the SL range are severely blended with surrounding PAH emission. The dotted lines show 
where the SL, SH, and LH data have been joined. Cosmic ray artifacts are indicated by " x". 



ent than the canonical (thermal equilibrium) value of 3, 
as in the presence of shocks (Wilgenbus et al. 2000), or 
by the superposition of multiple layers along the line of 
sight with different physical conditions. From this stand- 
point we must view our spatially-limited measurements 
as mean indicators of the physical conditions in both the 
PDR and outflow. 

The factor ~2.5 in Tox between the PDR and the out- 
flow seems quite reasonable, for line emission in the out- 
flow arising from coUisional (dissociative J-shock) excita- 
tion. We favor the J-shock mechanism over C-shocks in 
the outflow, based on the S(0) line intensity and strong 
[Si 11] 34.815 /zm emission. In the outflow the [Si 11] 
and [Fe ll] line fluxes are 1.06 (±0.15) xlO'^o and 2.36 
(±0.38) X 10-21 W cm-2, respectively. In the PDR, the 
[Si 11] line flux is 7.92 (±0.13) xlO'^i W cm^^. Inter- 
estingly, the measured line intensities of the S(0) line in 
the PDR and outflow are quite similar, which might be 
explained by a photodissociating component in the out- 
flow. As for velocities, Edwards & Snell (1983) measured 
a maximum CO expansion of 10 km scc^^, associated 
with the outflow. In J-shock models (HoUanbach & Mc- 
Kee 1989) for a faster jet-related component, H2 lines 



arc not a good diagnostics of shock velocities, as their 
intensities are almost flat as function of velocity at the 
densities we expect (~ 10^ cm"^). But we know from 
the absence of [Ne 11] and presence of H2 that the veloc- 
ity must be less than 50 km sec^i. At the densities and 
low Tex we measure, we estimate J-shock velocities in the 
range of 10 - 30 km sec"!. 

The II2 spectrum (and thus the Tex we infer) for the 
II2 ridge can be explained by the standard Draine & 
Bertoldi (1996) models for a stationary PDR, and are 
similar to conditions in the S 140 PDR (Timmcrmann 
et al. 1996). This does not rule out contribution from 
non-radiative processes that are reasonable in this en- 
vironment, where either or both of BD +65°1638 and 
BD +65° 1637 are sweeping material outwards by me- 
chanical wind momentum, to form some II2 by coUisional 
excitation. Indeed far infrared spectroscopic observa- 
tions of the LKHa 234 region with multiple pointings 
using the ISO-LWS (Lorenzetti et al. 1999) already in- 
dicate that the PDR is the dominant component for the 
II2 fluorescent emission, but the presence of C-shocked 
material cannot be ruled out. Our observations lack the 
spatial information with additional pointings to make a 
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comparison with a more complex, clumped PDR model 
(e.g., Burton et al. 1990) as done by Lorenzetti et al. 
(1999; see their Fig. 6). It is anyway useful to note 
that coUisional excitation, if present in a clumped envi- 
ronment, cannot be done in cooling zones behind disso- 
ciative J-shocks (cf. HoUenbach & McKee 1989), where 
Tex is lower than 300 - 400 K at our densities, but 
more likely in C-shocks with pre-shock H2 volume den- 
sity, magnetic field, ionization fraction, and shock speed 
conditions found, for example, in GGD37 in Cepheus A 
West (Wright et al. 1996). This region is also a site of 
intermediate mass star formation with a massive outflow, 
and bow-shock like clumps and filaments of H2 emission 
(Hartigan et al. 1996), exhibiting very comparable H2 
rotational line fluxes (and Tox — 700 K). Aperture aver- 
age column densities are one order of magnitude higher 
than in the NGC 7129 H2 ridge, which gives us a corre- 
spondingly lower H2 mass of ~ 1O~^M0. 

4. SUMMARY AND CONCLUSIONS 

With the IRS we have obtained the first mid-IR spec- 
troscopy of regions around LkHa 234, in the recogniz- 
able outflow and in a nearby ridge of II2 emission. The 
measured intensities of the PAH bands in the II2 ridge 
allowed us to infer the strength of the UV radiation field, 
which is comparable to that estimated by Boulanger et 
al. (1998) in NGC 7023. The rotational H2 emission line 
ratios infer an (LTE) aperture average excitation temper- 
ature of 657 K and column density of 1.2 xlO^^ cm~^, 
quite similar to the warm-Il2 regions of Cep A that are 
well matched by formation in non-dissociative C-shocks. 



The cooler II2 emission in the outflow is consistent with 
formation by coUisional excitation in J-shocks, with ve- 
locities of 10 - 30 km sec~^. A photodissociating compo- 
nent is suggested by similarity of S(G) line intensities in 
the outflow and PDR/II2 ridge. There is no significant 
contribution of thermal emission from warm dust in the 
outflow. 

The origin of the PDR cannot be localized with cer- 
tainty from the IRS observations presented here. Loren- 
zetti et al. (1999) presented ISO-LWS observations of 
[C 11] 158 /im, [O i] 63 /im, and [O i] 146 ^m emission 
whose spatial morphology and derived densities and UV 
field strength across the nebula favor BD-|-65° 1637, ar- 
guing that this is consistent with observations by Fuente 
et al. (1998) of a ^'^CO cavity created by the young star. 
Miskolczi et al. (2001) have subsequently shown with 
^^CO and ^'^CO maps that the cavity is instead centered 
on the older star BD+65° 1638, which Matthews et al. 
(2003) recently find in 21 cm maps to be surrounded by 
an extensive photodissociated H2 and a small H 11 region. 
On the other hand, our detection of warm material at 16 
^m around BD-|-65° 1638 clearly shows that this region 
is not free of material, and though BD-|-65° 1638 is not 
an embedded object, further spectroscopic observations 
are needed to ascertain its nature. 
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manuscript. 
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